Exploring correlated ID Bose gases from the superfluid to the Mott-insulator state by 

inelastic light scattering 
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We report the Bragg spectroscopy of interacting one-dimensional Bose gases loaded in optical 
lattices across the superfluid to Mott-insulator phase transition. Elementary excitations are created 
with a non-zero momentum and the response of the correlated ID gases is in the linear regime. The 
complexity of the strongly correlated quantum phases is directly displayed in the spectra which ex- 
hibit novel features. This work paves the way for a precise characterization of the state of correlated 
atomic phases in optical lattices. 
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Cold atomic gases loaded in optical lattices have been 
routinely used during the past few years as a versatile 
and powerful experimental system to study many physi- 
cal problems [ij. In particular, they allow to realize and 
manipulate strongly correlated quantum phases [3, [3, 01 
and constitute a promising candidate for implementing 
quantum information processing and quantum simula- 
tion schemes To achieve those goals, a corner stone 
consists in a precise characterization of the correlated 
gaseous phases. 

As in similar problems of condensed-matter, the pres- 
ence of strong correlations makes it hard to draw a com- 
plete picture, both from the experimental and theoreti- 
cal point of view. The implementation of spectroscopic 
probes, such as angle-resolved photo-emission spec- 
troscopy for high-Tc superconductors where electron- 
electron correlations play a major role [f| , are crucial. For 
the Mott-insulator phase created with trapped atomic 
gases, experiments have demonstrated the existence of 
a gap in the spectrum 0, 01 and have investigated the 
shell structure K , the spatial order via noise correla- 
tion techniques [0] and the suppression of compressibil- 
ity an experimental probe giving direct access 
to important information, such as the temperature or 
the elementary excitations on which the dynamical prop- 
erties of the many-body system depend, is still miss- 
ing. It has been proposed that inelastic light scattering 
(Bragg spectroscopy) performed at non-zero momentum 
in the linear response regime could provide such a tool 

0, El EE, 111. 



In this letter we report the measurement of the lin- 
ear response of interacting one-dimensional (ID) gases 
across the superfluid (SF) to Mott-insulator (MI) tran- 
sition. The ID gases are loaded in an optical lattice 
whose amplitude drives the atomic sample from a SF 
to an inhomogeneous MI state. The elementary excita- 
tions are created at a non-zero momentum using two- 
photon Bragg transition. The presence of an additional 
mode to the phonon mode in the SF state is suggested. 
In the inhomogeneous MI state, multiple resonances are 



observed which give information about the particle-hole 
excitation energy, the inhomogeneity of the trapped sys- 
tem and which exhibit novel features at low energies that 
could be related to the temperature of the atomic sam- 
ple. From the continuous modification of the spectra we 
also get quantitative information on the critical lattice 
amplitude for entering the MI phase. 
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FIG. 1: a Experimental sequence: after loading the 3D BEG 
in a 3D optical lattice {s± = 35, Sy) with an exponential ramp 
(time constant 30 ms), the ID gases are excited with the 
Bragg beams for a time Ate. The optical lattice is then lin- 
early ramped down to s± = Sy = 5 in 15 ms. After a waiting 
time of 5 ms, both the optical lattices and the magnetic trap 
are abruptly switched off and an absorption image is taken 
after a time-of-flight of 20 ms. b Typical absorption image 
in the absence of the Bragg excitation, c-h Spectra in the 
lowest energy band measured across the SF to MI transition. 
Lines are guides to the eye. Note the drop in the amplitude 
of the response (vertical scale). 

We produce an array of independent ID gases by adia- 
batically loading a three-dimensional (3D) Bose-Einstein 
condensate (BEC) of ^"^Rb atoms {N ~ 1.5 x 10^, chem- 
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ical potential /i3D — 740 Hz) in two orthogonal optical 
lattices (Al = 27r/gL = 830 nm) at a large amplitude 
V± = s^Eji = SSSr where Ej^ ~ /2mX\, h being 
the Planck constant and m the atomic mass. The ID 
gases can be considered independent since the resulting 
tunneling rate in the transverse plane is 0.75 Hz much 
smaller than the inverse time scale of the experiment. 
A crucial quantity indicating the regime in which a ID 
gas lies is the ratio of the mean-field interaction energy 
to the kinetic energy needed to correlate particles on a 
distance 1/niD, riio being the atomic density |15i]. This 
ratio is 7 = mgYDl?T?nm, giu being the inter-atomic 
coupling in the ID gas [l6[. For 7 -C 1, the ID gas is 
in the mean-field regime while for 7 3> 1 it enters the 
Tonks-Girardeau regime [13]. In our experiment 7 0.6 
[3] so the correlations in the ID gases are stronger than 
in the mean-field regime. To drive the ID gases from 
a SF state to a MI state, a longitudinal optical lattice 
(Al = 830 nm) with an amplitude Vy — SyEji is added 
along the axis of the ID gases. 

The spectroscopic study of the ID gases is performed 
using a two-photon Bragg excitation [19|. Once the ID 
gases are loaded into the optical lattices, two off-resonant 
laser beams (As = 780 nm, red-shifted by 250 GHz from 
the atomic transition), detuned from each other by a tun- 
able frequency difference are shone onto the atoms dur- 
ing a time At^ = 3 to 6 ms. The momentum transfer go 
by the Bragg excitation depends on the angle between the 
two beams while the energy transferred is given by the 
frequency difference v between them. The angle between 
the Bragg beams has been set up to transfer a momentum 
along the y direction smaller but close to the edge of the 
first Brillouin zone where the response in the Mott state 
is predicted to be the largest (lol . UH . The momentum 
transfer qq along y is kept constant to go = 0.96(3)gL 
[2o| . The spectra are measured by tuning v and moni- 
toring the response of the ID gases. The experimental 
sequence is described on FiglT^. From the interference 
density pattern after expansion (Fig.[Tj3), we fit the cen- 
tral peak with a Gaussian function from which we extract 
the RMS width a. We have checked experimentally that 
the increase of a induced by the Bragg beams is propor- 
tional to the energy transferred [lH . The response of the 
ID gases lies in the linear regime: on resonance (in the 
SF and in the MI regime) a increases linearly both with 
the Rabi frequency rie/STr — Vs/'^h and the duration 
AtB of the Bragg pulse, Vb being the amplitude of the 
moving lattice created by the Bragg beams; in addition 
Vb is at the one percent level of the amplitude Vy of the 
longitudinal lattice (for all Sy values, Vb < 0.15i?R). 

SHghtly different parameters (f^B, Aie) have been used 
across the SF-MI transition to obtain a good signal-to- 
noise ratio. In order to compare all the spectra, the mea- 
sured increase of a is rescaled with the parameter ile Ate 
and with the total atom number 2j] . This rescaled quan- 
tity is referred to as the amount of excitation and its 



plot versus the relative detuning 1/ as the excitation spec- 
trum. Our excitation scheme is different from previous 
techniques using lattice modulation [3, 0] for which the 
momentum transfer is zero and the modulation of the 
lattice amplitude is large (~ 25 %). 

SF-MI cross-over — The array of ID trapped gases 
undergoes a cross-over from a SF to an inhomogeneous 
MI state as the amplitude Sy of the longitudinal lattice 
increases In the absence of the lattice [sy = 0), the 
ID gases are superfluid at low enough temperature, a 
property we have checked measuring the frequency ratio 
of the breathing mode to the dipole mode 2J| . For large 
Sy [Sy > 10), the ID gases are deep in the inhomogeneous 
MI state 0]. This quantum phase transition appears on 
the spectra shown on Fig. [TJ for low Sy the response ex- 
hibits a single broad shape (see Fig.[Tt-d) while, for large 
Sy the spectrum shows a complex structure with multiple 
resonances (see Fig.[T^-h). In addition, the amplitude of 
the response drastically drops as the atomic system en- 
ters the insulating regime and goes deeper in the Mott- 
insulating phase (note the vertical scale on Fig.[T]F-h). 



a 




3500 


b 








^ 3000 










X 

^ 2500 










T3 
■5 

CO 2000 


M 

i 








1500 


, + - 

• 






1000 


I 





Amplitude s (Er) 



Amplitude s (Er) 



FIG. 2: a Resonant frequency in the lowest energy band. The 
solid line corresponds to the mean-field calculation in a single 
ID gas for the experimental trapping frequencies and atom 
number A^atom = 250. b RMS width of the resonance in the 
lowest energy band. The blue area is the expected position of 
the transition for trapped ID gases extracted from [25| and 
corresponds tol7/J = 8.5 — 9. 

The spectroscopy of the ID gases across the transition 
allows to mark the boundary between the SF state and 
the inhomogeneous MI state. In a periodic potential, the 
energy distribution of the states is captured by the Bloch 
bands picture. The energy bands of the lattice flatten as 
the amplitude of the periodic potential is increased. For 
a momentum transfer close to the edge of the first Bril- 
louin zone gi,, the flattening of the energy bands implies 
that the resonant frequency of the transition within the 
lowest energy band decreases as Sy increases. In addition, 
the tunneling J along the axis of the ID gases decreases 
with Sy implying a loss of phase coherence and a larger 
distribution of quasi-momenta populated. Approaching 
the phase transition, the spectrum width becomes equal 
to the energy width of the band and decreases with Sy 
as the band flattens. These two effects - the decrease 
of the resonant frequency and the width of the transition 
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within the lowest energy band - hold for the dephased ID 
gases before the phase transition. Once a Mott domain 
appears, new resonances are expected within the lowest 
energy band, enlarging the spectrum since the excitation 
frequencies of the MI state (~ U and ~ 2U, see below) 
are larger than but close to that of the SF state. We 
have fitted the lineshape of the spectrum corresponding 
to the transition within the lowest band with a single 
Gaussian function as the system approaches the phase 
transition. The resonant frequency and the width are 
plotted in Fig. [21 The width clearly exhibits a minimum 
which we attribute to the appearance of a MI domain. 
This position, in the range U/J = 8 — 10 (see Eq. [Ij, is 
in agreement with recent Monte-Carlo simulations pre- 
dicting J7/ J 8.5 — 9 for trapped ID gases [25l |. 

Correlated superfluid regime — In any periodic system, 
an excitation at a given quasi-momentum is associated to 
multiple resonances corresponding to the different energy 
bands. We restrict our discussion to the excitations to- 
wards the first two bands as shown on Fig. [3ti. In the 
absence of the lattice {sy = 0) and for the regime of inter- 
actions of the ID gases (7 ~ 0.6), the resonant frequency 
of the Lieb-Liniger model is indistinguishable from the 
mean-field solution within our experimental resolution 
[2^. Therefore the modification of the spectrum com- 
pared to the mean- field regime mainly consists in the 
appearance of a tail at low energies 27|. We compare 
the Bogoliubov band calculation in the periodic poten- 
tial with the measured resonant frequencies. We find a 
good agreement for the lowest energy band while the res- 
onant energy of the transition towards the second band 
is shifted towards higher energies as the atomic sample 
approaches the SF-MI transition (see Fig. [8)3)). 
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FIG. 3: a Spectrum of ID gases in the SF regime {sy — 
4). The solid line is a guide to the eye. The shaded areas 
correspond to the width of the two lowest bands as shown in 
the inset, b Resonance frequencies of the first two bands in 
the SF regime. The lines are mean-field calculations for a ID 
BEG with the experimental parameters (A'^atom ~ 250). 

As stated before the presence of correlations in the ID 
gases should manifest itself on the shape of the spectra 
rather than the peak resonant frequency. At Sj^ = the 
width we measure (2.5 kHz, see Fig. [!}:) is much larger 



than the experimental resolution (the spectrum width of 
the 3D BEC being 0.75(9) kHz for the same parameters 
of the Bragg excitation). We attribute this observation to 
the presence of correlations and thermal effects peculiar 



to the beyond mean-field ID gases [28|. For Sy > 0, the 



shape of the spectra in the SF phase becomes asymmet- 
ric towards high energies (see Fig. [3^) . For comparison 
we have performed the spectroscopy of a 3D BEC loaded 
in an optical lattice at Sy = 6 measuring a symmetric 
profile. Numerous theoretical works predict that in ad- 
dition to the usual phonon mode, an extra gapped mode 
is expected in strongly correlated superfluids whose am- 
jlitude is maximum at the edge of the Brillouin zone 
I3I 14, 2^. The momentum transfer qo = 0.96qL is thus 
appropriate for this observation. The presence of two 
excitation modes not resolved could support the asym- 
metry of the spectrum. On the contrary we can exclude 
other possible explanations. First, the asymmetric shape 
does not correspond to the algebraic tail of correlated ID 
Bose gases which is expected towards low energies [27| . 
Second, assuming that the longitudinal lattice dephases 
the ID gases implies a large quasi-momentum popula- 
tion that could support the asymmetric shape only for 
the transition towards the second band (see shaded areas 
on Fig. [3^). 

Inhomogeneous Mott regime — The response of the in- 
homogeneous MI phase exhibits a complex structure with 
multiple resonances (see Fig.[TJ-h) that we now describe. 
For large amplitudes Sy the ID gases can be described by 
the Bose- Hubbard Hamiltonian [301, 



H = -J 



<i,J> 



fi-C.) + ^ TLj [Uj 



(1) 

Here aj, are the creation and annihilation opera- 
tor of one boson at site i and rii — a\ai is the parti- 
cle number operator. The on-site interaction energy is 
given by U and the next neighbor hopping amplitude 
by J. ei is the energy offset experienced by an atom 
on site i due to the presence of the harmonic confin- 
ing potential. In the homogeneous MI state, the low- 
est excitation is a particle-hole (ph) excitation which 
has a minimum energy cost at zero momentum trans- 
fer AOj^ = U^l - 4(2no -f- l)J/U + (2 J/C/)2 ^ISi] with no 
being the filling factor (no = 3 at the center of our ID 
gases). Apj^, the so-called gap of the MI state, is smaller 
than U and asymptotically reaches U for large ratio U / J . 
We observe two peaks at frequencies larger than the gap 
Apjj. Since the excitations of a SF or a normal gas are 
expected below this threshold, we attribute these two 
peaks to the MI domains. The measured frequency of 
the lowest of these two peaks corresponds to the energy 
Aph((7o) for exciting a particle- hole at the finite momen- 
tum go E^nd it is plotted on Fig. [4^ as a function of the 
ratio U / J . The frequency of the second peak is twice the 
frequency of the first one with their ratio plotted in the 
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inset of Fig|4]a. Since the response of the atomic gases to 
the Bragg excitation lies in the Hnear regime, the peak at 
the energy 2Aph(go) does not come from non-Hnear phe- 
nomena. We attribute it to the inhomogeneity of the MI 
state due to the presence of the trap and a loading in the 
optical lattices which would not be fully adiabatic [31 1. 
Information about the inhomogeneity of the Mott state 
is thus directly accessible using Bragg spectroscopy. We 
now compare the particle-hole excitation Aph(go) mea- 
sured in the experiment to the case of a homogeneous MI 
at zero temperature T = (green solid line in Fig. 3^) 
Our measurements of Aph(go) are well below the 
prediction for the homogeneous case. This deviation can 
be related both to the inhomogeneity and the finite tem- 
perature in affecting the ph excitation energy for which 
no complete theoretical predictions exist so far. 
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FIG. 4: a Points correspond to Aph(go) measured in the inho- 
mogeneous Mott phase. The green solid line is Aph(go) calcu- 
lated for a homogeneous MI state at T = with no = 3, the 
dotted blue line is Apj^ and the dashed blue line is U. Inset: 
ratio between the two resonant energies in the Mott phase, 
b Low energy part of the spectrum measured at Sy = 13. 
The green solid line is a guide to the eye. The shaded area 
corresponds to 4 J — 12 J and the vertical dotted line is Api-j. 



Compared to the spectra using lattice modulation 7], 
novel excitations are observed in the MI state at low ener- 
gies for all the range U/J tested (see Fig.HJa). Since their 
frequency is smaller than A^^ they can not be attributed 
to an excitation from the ground state of the Mott state. 
Nevertheless, excitations are expected at these low fre- 
quencies in trapped systems at finite temperature. First, 
the excitation energy of the SF domains with a filling fac- 
tor smaller than no is ~ An^J for a momentum transfer 
close to the band edge [l3| ■ In our experiment we have a 
maximum filling factor tiq = 3 and we have plotted the 
energy range 4 J — 12 J as a shaded area on Fig. [3]b. Sec- 
ond, a non-vanishing temperature T allows the popula- 
tion of several ph excitations which should be observed in 
the energy transfer induced by Bragg spectroscopy [ll|. 
The excitations at energies larger than 12 J could be ex- 
plained by the finite temperature of the system. Yet how 
to extract a temperature from the position and/or width 
of these low-energy excitations is not clear and we hope 
that these experimental measurements will trigger a new 
interest in working out such a possible relation. 



In conclusion we have demonstrated the possibility to 
monitor, in the linear regime, the response of correlated 
quantum gases to an excitation at non-zero momentum. 
We have studied interacting ID gases from the SF state 
across to the inhomogeneous MI state. The complexity 
of such correlated quantum phases appears in the spec- 
tra from which important information is directly acces- 
sible. In the future, a systematic investigation of the 
observed novel features combined with new theoretical 
inputs should allow to adress open questions regarding 
those strongly correlated gaseous phases. 
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